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ABSTRACT 

In this work an extensive set of ground-based measurements of the diurnal 
variation of medium frequency radio wave absorption and virtual height is 
analysed in terms of our current understanding of the D- and lower E-region 
ion production and loss process. When this is done a gross discrepancy arises, 
the source of which is not known. 
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A STUDY OF ELECTPON DENSITY PROFIT.FS IN RELATION TO 
IONIZATION SOURCES AND GROiND-BASED 
RADIO WAVE ABSORPTION MEASUREMENTS, PART I 

1. INTRODUCTION 

In recent years computer modelling of photo-chemical processes has become 
one of the primary tools for advancing our uiKierstanding of the lower ionosphere. 
Since the early work of Keneshea (1967) computer models have been continually 
updated and expanded as more and more reaction rate coefficients and cross 
sections have become known from laboratory measurements. Unfortunately in 
many cases it is not possible to adequately test the predictions of a particular 
time-dependent chemical model by comparison with results from rocket-bome 
in situ measurements. This is because, from a practical standpoint, rocket 
measurements cannot be made frequently enougii to accurately trace the time 
dependence of the experimental quantity. However, as will be brought out in 
this report, our knowledge of at least one of the more important quantities, namely 
the electron density profile and its time dependence becomes much more precise 
when the results of in situ measurements can be combined with quantitative 
ground-based measurements of radio wave absorption. 

Extensive multi-frequency absorption data have been obtained for the equa- 
torial region over a long period of years by Gnanalingam (1972). Near the 
equator the anomalous day-to-day variability so characteristic of the mid-latitude 
winter mesosphere seems to be either absent or at least greatly reduced. It is 
well established that mesospheric pressures, densities and temperature undergo 
significant variation, both day-to-day as well as seasonal. The magnitude of 
both these variations however are smallest at lowest latitudes. This meteorological 
stability together with the presumed absence of energetic particle precipitation 
makes the equator the preferred location from which to study the build-up and 
decay of the lower ionosphere when it is under a simple solar zenith angle 
control. 

With the large body of absorption data available for the equator as our control, 
we propose two objectives: (1) to deduce from a consideration of ion-pair pro- 
duction and loss processes the diurnal variation of the D- and lower E-region 
electron density profile for the conditions of an undisturbed sun at solar cycle 
maximum: and (2) to deduce the noontime electron density profile for the con- 
ditions of a moderately disturbed sun. Objective (2) will be the subject of a 
separate report which we hope to publish later. In the present work we take up 
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objective (1) which we pursue as follows: We first construct *he best possible 
noontim'" reference electron density profile by combining the results of a rocket- 
borne exij: ment with the results of a series of ground-based multi -frequency 

absorption measurements. Next, from our current knowledge of ionization 
sources, wt calculate the noontime profile of the ion-production function. These 
two seeps immediately define the noontime profile of t the effective electron 
recombination coefficient. We now assume that this profile of is time 
independent. From the calculated solar zenith angle dependence of the ion- 
production function, we can then predict the decay of the electron density profile 
as the solar zenith angle moves toward sunset. Finally our predictions are 
tested against highly accurate measurements of the diurnal variation of radio- 
wave absorption. The work is organized as follows: - 

In section 2 we compile an atmospheric model adiich is representative of the 
equatorial region during equinoctial months. The model up to 95 km is taken 
directly from the tables given by Groves (1971) and extended to 160 km by means 
of the data given by Jacchia (1971). F the important minor consituent nitric 
oxide, we use the profile of Meira (1970). 

In section 3 we describe our electron collision frequency profile which is 
derived from the atmospheric model of section 2. Included in this collision 
frequency profile are first order correction terms for the effects of electron 
collisions with ions and with oxygen atoms. 

In section 4 we consider the solar radiations that govern the ion production 
function in the lower ionosphere. It is found that our knowledge concerning the 
variation of the X.U.V. flux with solar activity is rather poor. This situation is 
unfortunate since it makes it difficult to extrapolate flux measiurements to dif- 
ferent levels of sol ' activity. Nevertheless we proceed in section 5 to calculate 
our ion production function versus altitude and solar zenity angle using adiat we 
think are the most appropriate values of solar X.U.V. flux for the conditions of 
an undisturbed sun at solar cycle maximum. 

In section 6 results from a noontime rocket measurement at Thumba, India 
are combined with the results of multi -frequency absorption and virtual hei^t 
measurements made at Colombo, Ceylon on 72 selected days. This combination 
of in situ and ground-based measurements is used to construct a reference 
electron density orofile for a solar zenith angle y = 10“ for the conditions of an 
undisturbed sun at solar cycle maximum. By combining this reference electron 
density profile N^ (ref.) with t<.e ion production fvmctlon Q(10“) calculated for 
V - 10“, we can define an effective recombination coefficient a^ff=Q(ir“>/N^ (ref.) 
Our altitude profile of and its implications are discussed in section 7. 
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In section 8, starting with the reference electron density profile as an 
initial condition, the solar zenith angle dependence of the electron density N, is 
calculated by numerically solving the familiar equation: 

dN/dt =Q(v) - 

where Q(v) is the ion production function computed in section 5. 

From the electron density profiles of section 8 together with the collision 
freqxiency model of section 3 we can calculate the solar zenith angle dependence 
of radio wave absorption and virtual reflection height. This is done in section 9 
where we find that vdien the calculated values are compared against the measured 
diurnal variations, a gross discrepancy is revealed. Until this discrepancy is 
resolved, which we have not been able to do, our understanding of the basic 
processes of the lower ionosphere must be considered unsatisfactory. 

F'^ ally. in section 10, for the sake of completeness, we construct an ad hoc 
set of Ng (z.v) profiles which yield agreement with the measured diurnal varia- 
tion of absorption and virtual height. It is believed that this set of profiles can 
provide a fairly stringest test for any theory of the lower ionosphere. 


2. ATMOSPHERIC MODEL 

Althov^h we shall calculate the ion production function only to 110 km, we 
need an atmospheric model to at least 160 km to determine the attenuation of 
solar flux. The time independent atmospheric model which is used throughout 
the present work is based on the lower altitude data of Groves (1971) for 10” 
latitude at equinox, and extrapolated to 160 km by means of the Jacchia (1971) 
model for an exospheric temperature of 1100 degrees Kelvin. This value of 
exospheric temperature was chosen so as to be consistent with a 10.7 cm solar 
flux level of 135 units, i.e. the solar condition under which our radio wave absorp- 
tion rueasurements were obtained. Details of the model together with tabulated 
values are given in appendix A-1. 


3. COLLISION FREQUENCY MODEL 

The collision frequency model used in the present work is based on the 
recent review of Phelps (1972). For the ionospheric temperatures with which 
we are concerned, monoenergetic electron collision frequencies for the major 
constituents N^ and O^ can be combined conveniently into the expression 
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= (6.5 ±0.7) X 10®p„„ . 

where is the gas pressure in newton/m^. Above 90 km it is advisable to 
include first order correction terms for the effects of electron collisions with 
ions and electron collisions with atomic oxygen. The electron-ion collision 
frequency (Phelps, 1972) is given by 

^ie = 3.6N^T;*-5 x ln(2x 10^Te‘- 5/Ne°- 

where N, and are the electron density (om~^) and temperature C*K) respectively. 
To obtain the electron-atomic oxygen collision frequency , we use a theoretical 
expression given by Banks (1966) for the average momentum transfer cross section 


= (3.4 ± 1.0) X I0"‘«cm2 


For we then have 


= [0] X {1.88 X 10‘1°T^-®} (1 ± .3) 


where [O] is the density (cm~^) of atomic oxygen. 

The total collision frequency v = ^ to components are listed 

in appendix A-2. 

For the purpose of computing radio wave absorption, we use in the conductivity 
integrals the value of collision frequency given by = i^p + ^ . + Strictly 
speaking this is incorrect since it ignores the different energy depen^noe of the 
three components of v. However, in our case the radio waves are all reflected 
from altitudes below 105 km, where and are small correction terms. 
Hence the neglect of their different energy dependence should not lead to any 
significant error in our calculated values of absorption. 


4. SOLAR RADIATIONS 

A considerable volume of literature has been written concerning the solar 
radiations that are effective ionization sources for the D- and lower E-regions of 
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the ionosphere; see for example the reviews of Ohshio et al (1966), Ivanov- 
Kholodnyi and Nikol'skii (1969), Swider (1969), Aikin (1972). For the altitude 
region of concern in the present ivork, i.e. from 60-110 km, the major ionizing 
radiations are 

(a) the h 3 rdrogen Lyman-a line (1215. 7A). 

(b) the hydrogen Lyman-/3line (1025.7A) together with the nearby continuum 
in the 1025- 990 A interval, 

(c) the carbon m line at 977A and 

(d) soft X-rays in the wavelength interval loO-lA. 

The intensity of a particular radiation can be expected to vary considerably with 
solar activity. This is particularly true for wavelengths less than lOA. In 
appendix A-3 we gtve a orief survey of the present state of knowledge concerning 
these solar fluxes. Based on that survey the flux values used inthe present work 
are listed in table I for the conditions of an undisturbed sun at solar cycle 
maximum. 

In Fig. 1 the energy flux values for the range 100-lA are normalized to unit 
energy interval of leV ( o f/ in table I) and plotted versus energy. For the 
energy range corresponding to 100-30A (i.e. 124 < U < 413 eV) our tabulated 
fluxes are based primarily on the most recent values given by Manson (1972) 
which have a claimed accuracy of ±25%. Below lOA (i.e. U > 1240 eV) our tabu- 
lated flux values are obtained by fitting satellite measurements of the 1-8A 
energy flux to a bremsstrahlung t)rpe spectrum for a coronal temperature of 
3.5 X 10** degrees K. (This is discussed in appendix A-3.) For the range 30-10A 
(413 ^ U ^ 1240 eV) our tabulated fluxes are given by the arbitrary straight 
line of Fig. 1 which connects Manson's values above 30A with the bremstrahlung- 
type spectrum values below lOA. Fortunately our choice of flux values in the 
30- 10 A range does not strongly affect our calculated value of the ion production 
function, as we shall see in the following section. 
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Estimated Solar Flux* 
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5. ION PRODUCTION FUNCTION 


The total ion-production function Q(z. v) at an altitude z and solar zenith 
angle y is the sum of many separate production functions, each for a particular 
wavelength or wavulength band, i.e. 

Q(Z, y) = q^(Z. y) 

" \ 


where we define 


Fy (z. y): value of the attenuated solar flux at wavelength K 

N. (z): number density of the j constituent of the atmosphere. 

cr^ : cross section at wavelength K for photon absorption by a 

^ molecule of species j. 


7^ : photoionization yield, i.e. the number of electron-ion pairs 

' formed when one photon of wavelength k is absorbed by a 
molecule of species j. 


5.1 Absorption Cross Sect! d ~onizatlon Yields 

The subject of photon absorption cross sections and ionization yields In the 
wavelength range 100-1 A has been elucidated in the previously mentioned review 
by Swider (lO" ))*. Based on the forrnuis^ quoted in that review together with 
additional tab.ilatlons given by Ohshio ei d (1966), the values useu in the present 
work are listed in Table II for the atmospheric constituents O, O,, N^ and NO. 


Recent work by Henke and Elgin (1970) ho» »ub»fontioted the formulas quoted by Swider (1969) 
which dote bock to the year 1933. 
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TABLE II 

Absorption Cross Sections and Ionization Yields 


K 

(A) 

a(0) 

(cm2) 

(cm2) 

^(Nj) 

(cm2) 

>(0) 


^(Nj) 

1 

8.43E"23 

1.76E"22 

1.03E"22 

340 

340 

384 

2 

6.40E“22 

1.31E"21 

7.47E~22 

158 

158 

179 

3 

2,20E"21 

4.29E"21 

2.46E"21 

99.2 

99.2 

11? 

4 

5.02E“21 

9.93E’21 

5.73E"21 

69.7 

69.7 

78.9 

5 

l.OOE-20 

1.90E"20 

1.10E*20 

52.5 

52.5 

59.4 

6 

1.67E*20 

3.20E'20 

1.87E“20 

41.1 

41.1 

46.6 

7 

2.65E"20 

4.94E"20 

2.76E"20 

33.3 

33.3 

37.7 

8 

3.80E“20 

7.22E”20 

4.27E‘20 

27.4 

27.4 

3".0 

9 

5.40E-20 

1.01E~19 

5.70E~20 

25.0 

23.0 

26.0 

10 

6.90E-20 

1.34E"19 

8.05E“20 

19.2 

19.2 

21.8 

10-15 

1.18E-19 

2.70E~19 

1.50E~19 

16.7 

16.7 

18.9 

15-20 

2.68E"19 

7.00E"19 

3,60E“19 

13.5 

13.5 

15.2 

20-23 

3.50E"19 

7.00E“19 

3.60E“19 

11.8 

11.8 

13.3 

23-31 

4.50E~20 

9.00E“20 

1.00E‘18 

10.1 

10.1 

11.4 

31-103 

see footnote* 





977 

0 

? 98E"18 

8.20E“20 

0 

0.62 

0 

990-1025 

0 

1.75E"18 

7.40E"20 

0 

f».413 

9 

1025.7 

0 

1.52E“18 

1.00E~21 

0 

00 

0 

0 

1215.7 

0 

9.00E"21 

6.00E“23 

0 

0 

0 

1215.7 


o-(NO) = 2.40L~18 


7 (N<'-1 = 0.81 



For the importont wavelength iiitervol 103-3 lA we take; 
otOj) ^ 2cr(0) = 2.27 X 10"23 \2-5 cfn“ 

= 1.33 < 10"^^ X.2.5 cm2 
y(0.1 =7(0) =1.02 +273 V’ 

7(NJ =1.15 +306 
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5^ Number of Molecules Along Line of Sight 

In order to calculate the ion production function at an altitude z for a given 
solar zenith angle v, we need to know how the phot'^n flux at a particular wave- 
length ^ is attenuated along its path from the sun to the altitude z. This we can 
write as 

L(y)-® 

N.(L) dL}- 


where F^(®) is the photon flux at the wavelength X. incident on the top of the 
atmosjdiere (table I), and dL is an element of the line of sight from altitude z 
to the sun. Assuming that our atmospheric model of section 2 is spherically 
83 mametric, the calculation of S<'z. v), the number of molecules along the line 
of si^t to thr sun, defined by 


F^(z. y) = F,(®) exp 


I 


S(z. V) = 


•L(V)- 


N. (L) dL 


is a straight-forward exercise in plane geometry, the details of udiich are given 
in appendix A-4. In that same appendix we also tabulate for several solar zenith 
angles the Stz.y) values for the three atmospheric constituents N^, and O. 

We have now covered all the material necessary to enable us to compute the 
major components of the ion production function for the lower ionosphere. Before 
we discuss the results of such calculations however, we shall consider briefly 
some additional sources of ionization that are included in our calculation of the 
total Q -function. 
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5.3 Additional Ionization Sources 


In addition to the q\ components computed from the flux values of table I and 
the cross section data of table II our total Q -function includes contributions from: 


5.3.1 Scattered Radiation 


As is well known scattered Lyman- a and scattered Lyman- are major 
ionization sources for the lower ionosphere during nigattime. What may not be 
obvious is that these sources are not completely negligible during daytime, 
especially at hi^ solar zenith angles near twilight. To calculate the q due to 
scattered Lyman alpha and scattered L 3 nr.an beta we assume as a first approxi- 
mation that these sources are isotropically distributed over the u|^r hemisphere 
of the sky. The calculation of the scattered radiation is then only 

slightly more complicated than the calculation of q(y, z) for direct radiation, in 
that involves an additional integration over a local zenith angle ^ , i.e. 




^7T/2 

Jo 


q(6^, z) s in ^ 


For the intensity of daytime scattered Lyman alpha we use the twilight value of 
15 kilorayleigh (Meier and Mange, 1970). The scattered flux incident upon the 
top of the atmosphere is then approximately 2% of the value for direct solar 
Lyman alpha flux. Assumii^ this same ratio of scattered to direct flux, we take 
the intensity of the daytime scattered Lyman beta to be 150 rayleigh. 

5.3.2 Metastable O* ( ) 

2 ' 


To account for the ionization of the metastable oxygen state O* ) by 
solar radiation near 1108A, we have used the formula given by Paulsen et al 
(1971). 


q*Cz, v) = ND*(z) {0.549 X 10“’ exp(-2.406x 10"»S^^(z, y) 

+ 2.614 X 10“^ exp(- 8.508 X 10”^®S„ (z, y',} 

°2 

where ND*(z) is theO* (‘a ) number density and is the munber ofO^ 
molecules along the line of sight. This formula tak^s account of the attenuation 
of the 1108A flux by the minor constituent carbon dioxide. The authors of this 
formula estimate that there is an uncertainty of a factor of two on q* due to 
uncertainties both in the 1108A flux and in the cross sections that are involved. 
The O* profile given by Huffman et al (1971) is listed in appendix A-1, table A- 1-2 
Although this O* profile is based on measurements near 70° solar zenith a,igle, 
no significant increase in this profile is expected at lower solar zonith angles, 
on theoretical grounds. 
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5.J.3 Energetic Particles 


V, ? assume that near the magnetic equator precipitated energetic electrons 
can be .neglected as an ionization source for tbs O- and lower E- region. How- 
ever we 1o include in our total Q -function a contribution to galactic cosmic rays. 
To make this calculation we use an empirical formula (Mitra, 1968) adiich for 
magnetic latitudes less than 50*, fits the cosmic ray ionization rates measured 
by Van Allen (1952). 

^ (z) = ND(zV(1.5x lOr*® cos"-* L^) • {1 - 3 X 10"^(F - 70)} 

where ND is the total molecular number density, is magnetic latitude and F 
is the 10.7 cm solar flux in units of 10~^* W m“* Hz'* . 


5.4 Q- Function Resulto 

At a particular altitude the relative contribution to the total Q function made 
by a component is a function ot solar zenith angle. This is illustrated in 
Fig. 2 and Fig. 3 which show the q j components for solar zenith angles of 10* and 
60* respectively. Numei ical values of our results together with a calculated 
r.m.s. error on Q(z) are tabulated for four solar zenith angles in appendix A-5. 
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Figure 2. Calculated Q-function and components for x = 
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Figure 3. Calculated Q-function ond components for x = 60° 








6. REFERENCE ELECTRON DENSITY PROFILE 


The electron density profile has not yet been determined at equatorial 
latitudes from groimd-based methods such as the partial reflection or cross 
modulation experiment. The few profiles that are available for the equator 
have been obtained from rocket measurements at Thumba, India (Kane 1969, 

1972, Somayajulu et al 1971, Aikin et al 1972). The profile obtained at the 
smallest solar zenith angle for an undisturbed sun (i.e. 14°, Kane 1969) was 
adopted as the basis for the 60-82 km segment of our present reference profile. 
The uppermost portion of our reference profile was determined from empirical 
formulas which give as a function of solar activity and zenith angle the height 
and density at the peak of the E-region. In the altitude region between 82 km and 
the peak of the E-region at 105 km, our reference profile was adjusted by trial 
and error until calculated values of absorption and virtual height agreed with 
the extensive multi-frequency absorption and virtual teight measurements made 
at Colombo, Ceylon by Gnanalingam (1972). Details of the procedure for con- 
structing our reference profile are given in appendix A-6. The final profile tc 
emerge from that procedure is shown as the solid curve in Fig. 4. Also shown as 
the dashed curve in Fig. 4 is an alternate reference profile allowed by the confi- 
dence limits on the ground-based measurements. We shall have occasion to refer 
again to this alternate profile later in section 9. 

1 

h 

(km) 

110 

100 

90 

80 

70 

60 


Figure 4. Reference electron density profiles for x= Solid curve N^(ref) 
is deduced from mean values of multi frequency L, h ' meosurements (table ill), 
while dashed curve N^(ref) is the profile ol lowed by the upper limit of error on 
the L,h ' meosurements. 
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TABLE in 

Model Profile Tested Against 72 Days of Accurate Absorption Measurements 
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between and in the regfion below 85 km is to be expected on the basis 

that our a, fj contains a negative ion contribution, i.e. 

= (1 + ^)(aj + v«i) ^ an 

where S. is the ion-ion recombination coefficient and -V = N_/N, is the ratio 
of negative ion density to electron density. However since definitive values of 
and K(z) are not available, no inference concerning our Q-function or the 
controlling nitric oxide profile can be drawn from our large values of a^ff be- 
low 80 km. In the altitude region above 90 km the shape of our profile sug- 

gests that here our Q-function is deficient. Either there exists an additional 
ionization source which we have not taken into acccamt, or we have underesti- 
mated one or more of the solar flux components. We shall return to these pos- 
sibilities later in secti'^n 11. 


8. ELECTRON DENSITY VERSUS SOIJ^R ZENITH ANGLE 

To obtain the electron density profile as a function of solar zenith angle Y, 
we numerically solve the following equation: 


•N;(z, V) - Q(z. v) 


_d 


{Njz. V)} 


where 


\f( ^ Q(10'’)/N2(ref). 

For y ^ 75° we find that the dN/dt term is completely negligible, and that the 
electron density profil at zenith angle y is generated from our reference profile 
according to the expression 


v) = N^(rcf) 


/( 


Q(z- v) 

Q(z, 10°)' 


From this relationi,.iip it can be seen that the zenith angle dependence of the N^ 
profile will not be overly sensitive to the absolute values of the solar fluxes 
since what we are primarily concerned with here are ratios of Q, Calculated N^ 
profiles for V = 40°, 60° and 75° are plotted and discussed in appendix A-7. 

Implicit in our construction of the Np(z,y) profiles is the approximation 
the ratio y ^ N./N^ is time independent (i.e. not a function of solar zenith 
angle). For the altitudes above 85 km where most of the absorption takes place. 
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aiid for solar zeniih angles y < 75°, this approximation concerning negative iods 
should be quite adequate for purposes of calculating radio wave absorption. 


9. RADIO WAVE ABSORPTION: CALCULATED VERSUS MEASURED 

Having derived the solar zenith angle dependence of by the method described 
in section 8, we should now be able to predict the diurnal variation in the absorp- 
tion and the virtual height of reflection of a radio wave vertically incident on the 
lower ionosphere. In calculating these quantities by means of the generalized 
magnetoionic ray theory, the electron density and collision frequency profiles 
were logarithmically interpolated for altitude laminations of 100 meters. How- 
ever since the complex refractive index changes very rapidly in the region im- 
mediately below the height of reflection h^ , the last one kilometer of the calculation 
was done for laminations of 10 meter thickness. The height of reflection h^ is 
defined for the ordinary mode as the altitude at which the plasma frequency equals 
the frequency of the exploring radio wave. To the absorption calculated by ray 
theory, we add the so called phase integral correction which, near the magnetic 
equator is relatively small ^Thorpe 1971). At our exploring frequencies this 
combination gives a good approximation to a full wave solution. We estimate 
that the accuracy on our calculation of absorption is better than ±1.0 db. 

To calculate the virtual height we note that, for an ordinary wave propagating 
vert' Jly near the magnetic equator, the phase refractive index n and the group 
refractive index n' are related, to a good approximation, by the expression 
nn* = 1. Hence, from ray theory, the virtual height is given by the integral: 



To this value of virtual height a phase integral correction should be added to 
give a result closely equivalent to a full wave treatment. Instead however, for 
ease of computation, we make use of an alteraative procedure (Titheridge, 1967) 
where n(z) in the above integral for h* is evaluated with the collision frequency 
set equal to zero. We estimate that this procedure yields a virtual height which 
is accurate to better than ± 0.5 km. 

The measured diurnal variation of absorption and virtual height (Gnanalingam, 
1972) is given in table IV for the exploring frequencies 2.0 and 2.2 MHz. These 
results were compiled from measurements on 15 selected days in the equinoctial 
months of 1968, 1969 and 1970 when the solar conditions met our criteria for an 
undisturbed sun at solar cycle maximum, (see section 4, table I). 
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TABLE IV 

Measured Absorption and Virtual Height 


Frequency 

(MHz) 

j 

Local Time 
(hours) 

V 

(deg) 

L 

(db) 

6L* 

(db) 

h» 

(km) 

Sh'* 

(km) 

2.2 

12:25 

7.6 

51.8 

1.5 

HI 

1.1 


13:25 

20.3 

49.6 

1.6 

mm 

1.2 


14:25 

34.7 

47.5 


Wsm 

1.5 


15:25 

49.4 

39.0 


107.3 

1.1 


16:25 

64.3 

29.2 


110.3 

0.9 


17:25 

78.9 

20.1 

3.8 

114.0 

4.6 

2.0 

12:35 

10.1 

56.4 


101.0 

0.8 


13:35 

22.7 

54.2 

mm 

101.1 

1.3 


14:35 

37.2 

51.1 

1.7 

103.5 

1.3 


15:35 

51.9 

42.4 

2.3 

106.1 

1.4 


16:35 

66.8 

31.6 

2.7 

109.1 

1.2 


17:35 

81.4 

21.3 

3.8 

113.5 

3.9 


*95% confidence limit 


In Fig. 6 and Fig. 7 the calculated values of absorption and virtual height 
(shown as circles) are compared with the measured values of Table IV. A strik- 
ing disagreement between calculation and measurement is immediately obvious. 

In particular, the calculated virtual height does not increase fast enough as v 
increases from 10° to 60°. For y greater than 60°, the h' increases too steeply 
in a manner suggestive of group retardation. In the case of absorption, the 
variation with solar zenith angle seems to agree with measurement up to v = 60°. 
However, this must be considered fortuitous because the same Ne profiles do 
not yield the observed h' - y dependence. At v values greater than 60°, the 
peculiar variation of absorption is a direct consequence of the too rapid increase 
in the h* value. A consideration of these discrepancies indicates that at altitudes 
above 95 km our values do not decrease fast enough with solar zenith angle. 

In the above analysis with (z, y) profiles have been generated from what 
should be regarded as the most probable (ref) profile, because it is this 
profile which is most consistent with the mean values of the multi -frequency L, 
h' measurements made at y = 10°. It is clear however that we could have con- 
structed other Ng (ref) profiles which ai'e consistent with the confidence limits 
of the measurements. We therefore consider the possibility that by starting 
from a different N, (ref) profile we can generate a set of (z.y) values which 
yield better agreement with the measured diurnal variations. However, by trial 
and error, the best we could do along these lines is shown by the crosses of Fig. 6. 
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The corresponding electj jri density profile is that labelled N’(ref) in Fig, 4. 
previously referred to in section 6. In Fig. 6 it will be noted that even with the 
most favorable reference profile N'(ref) we can get marginal agreement between 
calculated and measured L, h' only out to \ - 40° on the basis of our present 
knowledge of the Q-functions. The apparent agreement in h' at v = 70° ~ 75° 
is fortuitous, being merely the intersection of two dissimilar variations, one 
calculated and the other measured. 

For the sake of completeness, however, it is of interest to find a set of 
electron density profiles which do yield agreement with the L, h*(y) measurements. 
This WQ do in the following section, even tnough we are unable to generate such 
profiles from our present understanding o the Q-functions. The calculated values 
of L, h'(/l losuli'ug froru these ad-hoc profiles are shown as the triangular 
points in Fig. 6 r- 

10. (z, y) MO OR THE LOWER EQUATORIAL IONOSPHERE 

AT EQUINOX 

We have already drawn actention to the fact that our present approach, based 
on our knowledge of the Q-fonctions, is only marginally adequate to generate 
profiles that are consistent with the measured diurnal variations up to Y = 40°. 

At this zenith angle the highest altitude which enters into our analysis is the true 
height of reflection on 2.2 MHz (i.e. 98.2 km). It would therefore seem that the 
Ng profiles need to be modified only at altitudes above 98.2 km. However, to 
get better agreement with the measurements we find that the modification should 
begin at about 95 km. Tt is interesting to note that this <s the level at which the 
slope of the profile begins to depart noticeably from the slope of the 
profile (Fig. 5). In modifying the profiles above 95 km we are constrained by 
the well-established relationships concerning the zenith angle dependence of foE 
and h„E (see Appendix A-G). The final profiles to emerge from these modifica- 
tions are tabulated in appendix A-8 and plotted in Fig. 8. 

A,- interesting feature predicvod by ouv equatorial N^(y) model is the rela- 
tively small y dependence of election density in the 85-90 km altitude interval. 

A test of this prediction, which could oe a fairly simple ground-based experiment, 
would be very helpful at this stage. Also shown in Fig. 8 by the dashed curves 
are the inadequate Ne profiles generated by the Q-fur.etions. Significant dif- 
ferences between the two sets of profiles are immediately evident. The question 
naturally arises whether these differences could tell us something about the 
shortcomings in our theory uf ion-production smd loss. 
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11. DISCUSSION AND CONCLUSION 


In attempting to answer the above question we note that (Fig- 5) is sub- 
stantially less than at all altitudes above 90 km. This would seem to indicate 
that there a component missing from our Q-function. Given the correct set 
of N, profiles (Fig. 8) and assuming above 90 km, we could deduce 

such a missii^ component. However, when we do this we find that the missing 
Q component decreases with solar zenith angle much faster than the decrease 
associated with a photon source. This must mean either that ionization sources 
other than solar photons and/or that loss processes other than simple recombi- 
nation are operative at the 100 km level in the equatorial ionosphere. 

Clearly the problem deserves consideration beyond what has been attempted 
in the present paper. Until it is resolved, our understanding of the basic processes 
of the lower ionosphere will continue to be in doubt. 
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APPENDIX A-1 
Atmoapheric Model 


The time independent atmospheric model which is used throughout the present 
work is based on the lower altitude data of Groves (1971) for 10** latitude at 
equinox, and extrapolated to 160 km by means of the Jacchia (1971) model for an 
exospheric temperature of 1100 degrees Kelvin. In column 2 of table A-1-1, M 
is the giam molecular weight taken from Jacchia. The bfmpe^'atures, pressures 
and mass densities from Groves are given in columns 3, 4 and 5. The number 
density MD is computed from the relation ND = (MD x 10'^) x (6.023E23) x M, 
and is listed in column 6. (6.023E23 is Avogadro's number). Below 90 km the 
molecular nitrogen density in column 7 is given by Nj * 0.7811 x ND, while the 
molecular oxygen density in column 8 is given by O, « 0.20955 x ND. At the 
altitudes 90 and 95 km the percentage of Nj and is given in the footnote at the 
bottom of the table. For the altitude region of 100 to 160 km the Nj, O, and O 
densities are taken from the Jacchia tables but enchanced by a factor of 1.1727. 
The factor of 1.1727 is necessary to normalize the Jacchia densities to the data 
of Groves in the overlsqjping altitude range of 95 to 110 km. 

In table A-1-2 we give the logarithmically interpolated values of our 
atmospheric model for the altitude region between 60 and 110 km. Included in 
table A-1-2 are models of two minor constituents, nitric oxide and metastable 
oxygon O* (*A ). The nitric oxide model is taken from the measurements of 
Meira (1970), while the values are those given by Huffman et al (1971). 


27 



TABLE A-l-x 


z 

(km) 

M 

(gram) 

T 

(deg) 

P 

(mks) 

MD 

(mks) 

ND 

(cgs) 

(cgs) 

— 

O 2 

(cgs) 

0 

(cgs) 


28.96 

221 

250(1) 

394(-4) 

819(15) 

640(15) 

172(15) 



28.96 

233 

119(1) 

l78(-4) 


289(15) 

776(14) 


35 

28.96 

249 


823(>5) 

171(15) 

134(15) 

359(14) 


40 

28.96 


306(0) 

mm 

844(14) 

660(14) 

177(14) 


45 

28.96 

272 

162(0) 


433(14) 

338(14) 

907(13) 


50 

28.96 

274 

883(-l) 

112(-5) 

233(14) 

182(14) 

488(13) 


55 

28.96 

265 

473(-l) 

622(-6) 

129(14) 

101(14) 

271(13) 


60 

28.96 

246 

247(-l) 


728(13) 

569(13) 

153(13) 


65 

28.96 

227 

121(-1) 

186(-6) 

387(13) 

302(13) 

811(12) 


70 

28.96 

210 

568('2) 

942(-7) 


153(13) 

411(12) 


75 

28.96 

200 


435(-7) 


707(12) 

190(12) 


80 

28.96 

199 



395(12) 

309(12) 

828(11) 


85 

28.95 

200 

471(-3) 

821(>8) 

171(12) 

133(12) 

358(11) 


90* 

2L.83 

197 


364(-8) 


597(11) 

157(11) 

779(9) 

95** 

28J23 

193 


157(-8) 

335(11) 

258(11) 

606(10) 

169(10) 

100 


201 

386(-4) 



105(11) 

217(10) 

128(10) 

105 


229 

182(-4) 



433(10) 

803(9) 

742(9) 

110 


278 

980(-5) 



188(10) 

315(9) 

434(9) 

115 






879(9) 

136(9) 

262(9) 

120 






449(9) 

644(8) 

166(9) 

125 






250(9) 

338(8) 

111(9) 

130 






151(9) 

193(8) 

786(8) 

135 






970(8) 

119(8) 

581(8) 

140 






657(8) 

771(7) 

446(8) 

145 






464(8) 

524(7) ' 

353(8) 

150 






338(8) 

369(7) 

286(8) 

155 






254(8) 

268(7) 

236(8) 

160 






194(8) 

199(7) 

199(8) 


*ot 2 = 90km Nj =0.7848 xND; Oj = 0.2059 xND 
**ot 2 =95 km Nj =0.7694 x ND; Oj =0.1809 x ND 


28 













TABLE A- x-2 


Atm jsphevic Model 


z 

(km 

O 

(cm 


^'2 

<cm :) 

SO 

(cm 3) 

0*2 

(cm 3) 

6.000E1 



5.686E15 

2.600E8 

I.658E10 

6.100E1 



5.010E15 

2.2S0E8 i 

1 531E10 

6.200E1 



4.413E15 

2.000E8 

1.414E10 

6.300E1 


1.044L . 

3.891Eir> 

1.730E8 

1.107E10 

6.400E1 


9.11 iEl 

3.429E15 

1.520E8 

1.207E10 

6.500E1 


106EU 

3.022E15 

1 .320E8 

1-115E10 

6.600E1 


7.v»75E14 

2.637E15 

1.150E8 

1.030E10 

6.7O0E1 


6.17=E14 

2.302E15 

1.050E8 

9.514E9 

6.800E1 


5.3s .£14 

2.009E15 

8.HOOE7 

8 7 39E9 

6.900E1 


4.704E14 

1.753E15 

7.700E7 

8.i:9E9 

7.000E1 


4.1'»5E14 

1.530E15 

6.800E7 

7.500E9 

7.100E1 


3.5, jEI4 

i.3nEi:> 

5.850E7 

6.928E9 

7.200E1 


3.0HE14 

1.123E15 

5.100E7 

6.400E9 

7.300E1 



9.626E14 

4.450E7 

5.879i;9 

7.400E1 


2.2P£14 

8.248E14 

3.900E7 

5.400E9 

7.i>00El 


1 596E14 

7.067E14 

3.370E7 

5.144 ,5 J 

7.600E1 


1.606E14 

5.988E14 

2.930E7 

4.'‘uOE9 1 

7.700E1 


1.361E14 

5.074E14 

2.600E7 

4.59UE9 

7.800E1 


1.153E14 

4.299E14 

2.280E7 

4.300E-: 

7.900E1 


9,772E13 

3.643E14 

2.030E7 

4.041E9 

8.000E1 


8.280EI3 

3.087E24 

1.800E7 

3.800E9 

8.100E1 


7.002E13 

2.610E14 

1.640E7 

3.541E9 

8.200E1 


5.920E13 

2.207E14 

1.500E7 

3.300E9 

8.300E1 


5.006E13 

1 o6E14 

I.400E7 

3.040E9 

8.400E1 1 


4.233E13 

1.-J78E14 

1 1.340E7 

2.800E9 

8.500E1 1 

i 1 

i 3.579E13 

1.334E14 

1.350E7 

2.565E9 

8.600E1 

j 

3.034E13 

i 1.136E14 

K410E7 

2.350E9 

8.700E1 


2.571E13 

1 9.671E13 

1.550E7 

2.113E9 

8.800E1 


2.179E13 

8.233E13 

1.800E7 

1.900E9 

8.900E1 


1.847E13 

7.010E13 

2.250E7 

1.688E9 

9.000E1 

7.793E11 

1.566E13 

5.968E13 

2.«30E’? 

1.500E9 

9.100E1 i 

9.101E11 

1.295E13 

5.045E13 

j 3.320E7 

1.342E9 

9.200E1 

1,063EI2 

1.071E13 

4.265E13 

1 3.820E7 

1 .200E9 

9.300E1 

1.241E12 

8.858E12 

3.606E13 

4.330E7 

1.028E9 

9.400E1 

1.449E12 

7.327E12 

3.049E13 

4.900E7 

8.800E8 

9.500E1 

1.693E12 

6.060E12 

2.577EJ3 

5.450E7 

7.563E8 

9.600E1 

1.601E12 

4.934E12 

2.152E13 

6.050E7 

6.500E8 

9.700E1 

1.514E12 

4.018E12 

1..98E13 

6.700E7 

5.586E8 

9.800E1 

1.432E12 

3.272E12 

j 1.501E13 

7.400E7 

4.800E8 

9.900E1 

i.364El2 

1 2.664E12 

' 1.254E13 

8.100E7 

1 4.040E8 

l.OOOE^ 

I.281E12 

j 2.169EI2 

I.047E13 

9.000E7 

j 3.400E8 

1.010E2 

1.148E12 

! 1.778E12 

' 8.779E12 

9„500E7 

2.915E8 

1.020E2 

1.029E12 

1 1.457E12 

7.358E12 

1.000E8 

2.500E8 

X.030E2 

9.228E11 

1.195E12 

; 6.167E12 

1.040E8 

2.179E8 

1.040E2 

8.273E11 

9.793E11 

1 5.169E12 

1.090E8 

1.900E8 

1,050E2 

7.416E11 

8.028E11 

; 4 333E12 

1.120EH 

1.631E8 

1.060E2 

6.662E11 

6.659E1I 

3.666E12 

1.150E8 

' 1.400E6 

1,070E2 

b.984Ell 

5.523E11 

3.102E12 

1.180E8 

1 1.183E8 

1,080E2 

5.375E11 

4.f)81En 

2.625E12 

j 1.180E8 

1 1.000E8 

1.090E2 

4.828E11 

3.800E11 

2.221E12 

1.130E8 

j 9.055E7 

1,100E2 

4.337E11 

3.152E11 

1 1.879E12 1 1.M0E8 

1 8.200E7 
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APPENDIX A-2 

Electron Collision Frequency Model 


Throughout the present work the electron collision frequency is assumed 
to be independent of solar zenith angle. The total collision frequency v = v + 
^eo components are listed below, is the component due to 

electron collisions with and x'lj molecules, is the component due to 
'’oulomb collisions, and is the component due to electron- atomic oxygen 
collisions. 



blank NOT F1UC3) 
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TABLE A-2-1 
Collision Frequency 


^ ! 

(km) I 

p 

(sec *) 

(sec ') 

(sec *) 

1 

{sec *) 

60 

1.606E7 



1.6U6E7 

61 

1.392ET 



1.392E7 

62 1 

1.207E7 



1.207E7 

63 

1.046E7 



1.046E7 

64 

9.072E6 



9.072E6 

65 

7.865E6 



7.865E6 

r6 

6.761E6 



6.761E6 

67 

5.812B6 



5.812E6 

68 

4.996E6 



4.996E6 

69 

4.295E6 



4 295r6 

70 

3.692E6 



3.69.ii. . 

71 

3.133E6 



3.133E6 

7? 

2.6r9E6 



2.659E6 

73 

2.256E6 



2.256E6 

74 

1.915T:6 



1.915E6 

75 

1.625E6 



1.625E6 

76 

1.376E6 



1.376E6 

77 

1.166E6 



1.166E6 

78 

9.S75E5 



9.875E5 

79 

8.365E5 



8.365E5 

80 

7.085E5 



7.085E5 

81 

5.990E5 



5.990E5 

82 

5.065E5 



5,065E5 

83 

4.282E5 



4.282E5 

84 

3.621E5 



3.621E5 

85 

3.062E5 



, 3.062E5 

86 

2.595E5 

1 

i 

1 2.595E5 

87 

2.199E5 



2.199E5 

88 1 

1.864E5 



1.864E5 

89 

U580E5 



1.580E5 

90 

U339E5 

2.060E3 

1.469E2 

E361E5 

91 

1.130E5 

2.400E3 

1.907E2 

1 1.155E5 

92 

* 9.529E4 

2.796E3 

2.476E2 

1 9.833E4 

93 

1 8.038E4 

3.257E3 

3.215E2 

8.396E4 

94 

1 6.781E4 

3.794E3 

4.173E2 

7.202E4 

95 

1 5.720E4 

4.420E3 

5.418E2 

6.216E4 

96 

j 4.851E4 

4.197E3 

6.718E2 

5.338E4 

97 

4.114E4 

3.984E3 

8.329E2 

4.595E4 

98 

3.489E4 

3,783E3 

’ 1.033E3 

3.970E4 

99 

1 2.959E4 

3.592E3 

1.280E3 

3 446E4 

100 

2.509E4 

3.410E3 

1.587E3 

3.009E4 

101 

2.159E4 

3.098E3 

1.698E3 

2.63SE4 

102 

1.857E4 

2.814E3 

1.81HE3 

2 320E4 

103 

j 1.598E4 

2.557E3 

1.942E3 

2.048E4 

104 

1 1.375E4 

1 2.323E3 

2.076E3 

1.815E4 

105 

1 1.183E4 

1 2.110E3 

2.221E3 

1.316E4 

106 

1 1.045E4 

1.930E3 

2.105E3 

1.449K4 

107 

9.235E3 

- 1.765E3 

1.996E3 

1.30CE4 

108 

8.160E3 

1.614E3 

I.892E3 

1.167E4 

109 

7.210E3 

1.476E3 

i,794E3 

1.048E4 

no 

6.370E3 

1.350E3 

j K700E3 

9.420E3 
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APPENDIX A-3 
Solar Badiations 


HYDROGEN LYMAN- a (1215.7A) AND LYMAN- /?(1205.7A) 

Intensity measurements of the solar Lyman- a and Lyman- /3 emission lines 
made at different levels of solar activity are reviewed in Tables A-2.1 and A-2.2 
respectively. Also included in the tables are the average 10.7 cm solar flux 
values for the dates of these measurements. 

For the present work we need to know the mean intensities of these lines at 
a typical solar cycle maximum when the 10.7 cm fli x lies in the range 130-160 x 
10”^^ Wm”^ Hz“* (flux unit). In the case of the Lyman-a line, there are only a 
few reliable measurements reported for these solar maximum conditions. How- 
ever, several measurements are available for the solar minimum period 1964/65 
(Weeks 1967). In order to derive from the latter measurements a representative 
value of the intensity at solar maximum we examine the solar-cycle variation of 
both the Lyman- a and Lyman- /3 intensities. 

Inspection of the tables will show that, while there is a general increase in 
the intensities of these lines with solar activity, the scatter in the measurements, 
especially of the Lyman- ^ line, makes it difficult to evaluate the magnitude of 
this increase. A study of the solar-cycle variation of the sum of the intensities 
of several chromospheric emission lines, including Lyman-/3, has indicated an 
increase of the EUV flux by a factor of 1.5 for an Increase in the 10.7 cm flux 
by a factor 2, from 75 to 150 units (Hall et al 1969). These intensities were 
derived from rocket-borne measurements carried out during the period from 
August 1961 to November 1968. A similar relationship has been noted in the 
measurements made by a spectrometer aboard the satellite OSO 3 during a 
solar rotational period of 27 days In May 1967 when the opposite hemispheres 
of the sun exhibited very different levels of solar activity. The intensities of the 
Lyman- a and Lyman- /3 lines were observed to increase by factors of 1.3 and 
1.4 respectively when the 10.7 cm flux increased from 111 units to 201 units 
(Hall and Hinteregger 1970). These factors have an accuracy of approximately 
±0.1. Having regard to all the available observations, we estimate the solar- 
cycle factor to be 1.4 for both these lines. 

Applying this factor to the aversige value of the Lyman- a intensity at solar 
cycle minimum, which we estimate to be 3.9 ergs cm“^ sec"* from Table A-2.1 
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(after excluding extreme values), we obtain a figure of 5.4 ergs cm”^ sec”‘ or 
3.3 X 10 ph cm sec for the intensity at a typical solar cycle maximum. 

This figure is rather less than the value of 6.1 ± 0.45 deduced by Weeks (1967) 
for solar maximum, because his average value is weighted by the measurements 
made during the period of unusual sclar activity in 1957/58. 

The mean intensity of the Lyman- /3 line may be deduced directly from 
Table A-2.2 wl..ch contains an adequate number of measurements corresponding 
to this range of 10.7 cm flux. It is found to be 3.6 x lO^ ph cm”^ sec"* or 
7.0 X 10"^ erg cm“^ sec"' . 

THE CONTINUUM 1025-990 A 

Hinteregger et al (1965) report a value oi 2.7 x lo’ cm"^ sec”* for the 
photon flux in this wavelength band at solar minimuir conditions. By normalizing 
this value to the adjacent 1025. 7A Lyman- B flux also reported by Hinteregger et al 
for solar minimum, we obtain 4.2 x lo’ cm”^ sec”* for the value of photon flux 
that is used in the present paper. 


THE CARBON III LINE AT 977 A 

Published values of the intensity of this solar emission line measured by 
rocket- and satellite-experiments are to be found in Hall et al (1969) and Hall 
and Hinteregger (1970). For an undisturbed sim at solar cycle maximum we 
adopt the value of 4.4 x lo’ ph cm”^ sec”* (9.0 x io~* erg cm”* sec”*) recom- 
mended by Hinteregger (1970) in his review of the EUV solar spectrum. 


X-RAYS IN THE WAVELENGTH RANGE 100-lA 

For aeronomical calculations it is informative to divide this wavelength 
range into several X-ray bands which deposit their energy at different altitudes 
in the lower ionosphere. 

100-40A and 30-10A 

These bands are characterized by the presence of a number of spectral lines 
which have been investigated in some detail because of their importance in the 
study of solar physics. The greater part of the ionization produced by these 
bands occurs in the altitude range 95-120 Km, when the solar zenith angle is 
small. 
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The flux in the ba’*'*" 8-20A and 44-60A has been systematically monitored 
for a number of years by ion chambers on board the NRL SOLRAD satellites. 

Of the two bands, both of which produce a peak of ionization around 103 km 
{y = 10®), the 44-60A band is the vnore important from the point of view of 
E-region formation. 

The published flux data in this band have generally been derived from the 
ion-chamber currents on the assumption that the incident solar radiation may be 
represented approximately by a gray-body distribution at a temperature of 
0.5 X 10® °K (Kreplin 1961). Since the actual spectrum is marked by several 
lines, the published data are subject to some error. Even more serious is the 
fact that the 44-SOA photometers employed in the SOLRAD satellites respond 
with high efficiency to X-radiation of wavelength less than 20A (Kreplin and 
Horan 1969; NRL Report 6800, p. 39). Although the contribution to the photometer 
output from this radiation is relatively small when the sun is quiet, it is appreciable 
St solar maximum and frequently dominates the output during solar flares (Kreplin 
and Gregory 1966). The true values of the flux at solar maximum are probably 
much less than the uncorrected published data, which lie in the range 0.1 -0.6 
erg cm’^ sec’* for a non-flaring sun. 

Spectrometric measurements of the line intensities and continuum in this 
band (Manson 1967, 1968: Hinteregger 1970) have yielded integrated flux values 
that are significantly less than the SOLRAD data. Typical fluxes are 0.021 and 
0.030 erg cm’^ sec’* for the 40-50A and 50-60A bands respectively. In the 
present paper we have adopted the measurements given in a recent review by 
Manson (1972). Any systematic error in the calibration of the instruments 
from which these measurements were obtained has been estimated by him to be 
within ±25%. 


40-30A 

Solar radiation in this bond is the principal source of daytime ionization at 
altitudes around 90 km. Since the major part of the noon absorption of radio 
waves of frequency 1.5 - 2.5 MHz occurs at these altitudes, the 30-40A band of the 
solar spectrum is of particular importance m the interpretation of absorption 
measurements. Unfortunately the flux in the band has not yet been systematically 
monitored by detectors aboard satellites. 

Our present knowledge of this flux is mainly based on a few measurements 
made with rocket-gorne instruments, of the intensities of individual spectral 
lines within the band. The most prominent of the&e lines is the CVl line at 
33, 7A. Published measurements of the intensity of this line are shown in 
Table A-2.3. It will be noted that there is an inconsistency in the measurements. 
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For the same level of solar activity, as indicated by the 10.7 cm flux, the valuds 
differ by a wide margin, in some instances by nearly fa order of msignitude. 

The reasons for these discrepancies have not been resolved. In a search 
for an explanation Manson (1972) has examined the solar conditions that prevailed 
during each of several measurements carried out by the AFCRL group, the NPL 
group and the Culham group (Freeman and Jones 1970). His conclusion is that 
each high value of intensity reported and each remarkably strong image on a 
photographic spectrum which has been published for this line can be explained 
by flaring or 'near-Flaring' conditions. This line wr.3 only marginally detected 
when the solar spectrum was sampled under non-flai*ing conditions on 8 August 
1967 (10,7 cm flux = 143 units): and the background level too high to allow a 
meaningful integrated flux to be reported for the 30-40A band (Manson 1972). 

We have therefore adopted a value that is slightly higher than the flux reported 
by Manson for solar minimum. This value is also consistent with the flux in 
the neighbouring bands (see Figure 1). 


Wavelei^hs Less than 10 a 

X-rays of wavelength less than lOA produce a peak of ionization at an altitude 
of about 95 km (X' = 10°) when the sun is quiet, and at progressively lower alti- 
tudes as the hardness of the X-ray spectrum increases with solar activity. 

The flux in several bands in this short wavelength 'tali' ot the solar specirum 
has been systematically monitored for several years with instruments on board 
satellites. The 1-lOA flux values adopted by us are based on the SOLARD 9 
data. The reduced data are published In the solar Geophysical Data Bulletins 
as hourly averages of the flux in the bai«ds 1-8A and 8-20A. They are idso 
available in the form of minute-by-minute listings for these bands as well as the 
0.5 - 3A band, upon request to the NBL. A noteworthy feature of the measure- 
ments is that there is provision of the effects of particle interference on the 
ion chambers to be detected and filtered out of the measurements in the data 
reduction process (Kreplin and Horan 1969). 

In applying the data to aeronomical calculations it is Important to remember 
that the fliixes are derived by xising a gray-body spectrum. The assumed 
temperature is 2 x lO® K for the 1-8A detector and 10^ K for the 0.5-3A de- 
tector. Horan (1970) has re-evaluated the detector- current to energy-flxuc con- 
version constants for these detectors on the basis of a bremsstrahlung-type 
emission spectrum at different temperatures. The values obtained by him for 
the 1-8A detector are consistently less than the conversion constant derived 
from a gray-body spectrum at 2 x 10® K. Bremsstrahlung-type distributions 
at 3.5 X 10® K and 7.0 x lO® K give constants that are respectively 18% and a 
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factor of two less than the assumed value. Similarly, the constant for the 0.5-3A 
detector is 3-13% less than the assumed value, for temperatures in the range 
7-20 X 10® K. Horan (1970) has also shown how the temperature of the X-ray 
emitting regi^' ^.s, assumed to be isothermal, may be deduced from the ratio of 
the fluxes in the 1-8A and 0.5-3A bands. 

In the present work we have been unable to derive this temperature for an 
xindisturbed sun, using I'he method suggested by Ho an, because measurements 
of the 0.5-3A flux are not available for these solar conditions. The flux in this 
band from a non-flaring sun is generally below the SOLBAD 9 detector threshold 
even at sunspot maximum. We have therefore had to rely on estimates of the 
temperature obtained from other experiments (Neupert 1969, Pounds 1970) and 
have adopted a temperature of 3.5 x lO® K for the undisturbed sun at solar cycle 
maximum. The flux values have been calculated at lA intervals by fitting a 
bremsstrahlung-type spectrum at this temp>erature to the SOLBAD 9 measure- 
ment of 4 X 10"^ erg cm'^ sec~^ for the average 1-8A energy flux imder these 
solar conditions. 

For the moderately disturbed sun to be considered in a later Part n of this 
work the 1-8A flux lies in the range 2. 0-2. 5 x 10"^ erg.cm'^ sec'* . Even at 
this level of solar activity the 0.5-3A flux is barely metusursible by the SOLBAD 9 
detectors. However, by examining a large number of minute-by-minute listings 
of the fluxes in the 0.5- 3A and 1-8A bands, we have estimated the temperature 
of the X-ray active regions to be 7.0 x 10® K. 
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Table A-3.1 

Solar Lyman- a Measurements 
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Table A -3.2 

SoJar L 3 Tnan-/i Measurements 
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JRBfCEDING PAGi; NOT FZLMEZI 


APPENDIX A-4 

Number of Moleculr ■ Along Line of Sight 

In Tables A-4 we list for seven I solar zenith angles the number of 0» « 

and Nj molecules along the line of si^t to the sun (i.e. S(z»x) section 5.2). 
To calculate this quantity we assume our model atmosphere to be spherically 
symmetric. At an altitude z we construct 50 concentric shells, each of thick- 
ness r s 1 km. From the geometry of Fig. A-4 we have 



Figure A-4. Geometry for calculating S(z, y), the number of molecules 
olong the line of sight to the sun. 
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rJ = *^ + - 2Rj^3 cos(7t - Vj) (1) 

2R„-1^„ COS(7T - y„.j) (2) 

* 2tVl ~ Vn-l>^ - (*^n " *n-l) = ^ 


^n = R„_j cos(7T - v„_j) + ytVi - Y„ i)] ^ + (Rf- 


(4) 


sin sin(7T - y^) siny^ 

*. " K. 


(5» 


s in V sin Y 



*n ■ " 

I\) = ^ + z 

Rg = 6378.0 (radius of Earth at equator) 


( 6 ) 


( 7 ) 

( 8 ) 

(9) 

( 10 ) 



where 



is averse density In nth shell 


y = Z + n - 1 


(li) 


(12) 
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Table A-4.1 


Alt 


X = 40» 







6.00F1 

2.18718 

1.13721 

4.23721 

2.7<»;?ip 

1.4672 A 

5.44^21 

e.ioFt 

2.22^18 

88^20 

3.69^21 

2.85«P 

1.27F21 

4.74r21 

6.2071 

2.26K1R 

8.60720 

3.21721 

2.90718 

1.10721 

4.12721 

6.3071 

2.30718 

7.47F20 

2.79F21 

2.94F18 

9.59M0 

3.58JT21 

6.40^1 

2.33^18 

6.47W0 

2.42F21 

2.98718 

8.31720 

3.10IT21 

6.50^1 

2.35^18 

5.59W0 

2,09^21 

3.02^18 

7.18M0 

2.68K21 

6.60F1 

2.38718 

4.82^20 

1.80721 

3.Ct718 

6.19F20 

2.32JT21 

6.7071 

2.40718 

4.15720 

1,55F21 

3.08718 

5.33720 

1.99F21 

6.8071 

2.42718 

3.56720 

1.33721 

3.11718 

4.58720 

1.71721 

6.9071 

2.44718 

3.05720 

1.14721 

3.13718 

3.92720 

1.47721 

7.00^1 

2.46^18 

2.60720 

9*75^20 

3.16F18 

3.34720 

1.25721 

7.1071 

2.48718 

'2.21F20 

8, 31^20 

3.18718 

2.85^20 

1.07^21 

7.20F1 

2.49718 

1.88J720 

7.07MO 

3.20718 

2.42720 

9.09F20 

7.3071 

2.51718 

1.60720 

6.01720 

3.21718 

2.05720 

7.73M0 

7.40K1 

2,52^13 

1.36if20 

5.11720 

3.23718 

1.74720 

6.56F20 

7.5071 

2.53718 

i.isno 

4.33JT20 

3.25718 

1.47720 

s.seno 

7.60M 

2.54718 

9.70JP19 

3.67F20 

3.26718 

1.25720 

4.71F20 

7.70F1 

2.55iri8 

8.19P19 

3.10720 

3.27718 

1.0SE70 

3.99F20 

7.8071 

2.56718 

6.91719 

2.63723 

3.29F18 

8.88719 

3.38720 

7,90i?l 

2.57^18 

5.83719 

2.23^20 

3.30718 

7.49719 

2.86720 

8.C071 

2.58718 

4.92719 

1.88720 

3.31718 

6.32719 

2.42720 

8.1071 

2.59718 

4.14719 

1.60720 

3.32718 

5.32F19 

2*05;?20 

8.2071 

2.59718 

3.48F19 

1.35720 

3-.33i?18 

4.48719 

1.74720 

8.3071 

2.60i?18 

2.93J?19 

1.14720 

3.34F18 

3.76719 

1.47720 

8.4071 

2.61718 

2.46ri9 

9,70M9 

3.35iri8 

3.16719 

1.25720 

8.5071 

2.61^18 

2.06719 

8.22719 

3.35718 

2.65719 

1.06720 

8.6071 

2.62F18 

1.73719 

6.97719 

3.36718 

2.22719 

8.95719 

8.7071 

2,62«8 

1.44719 

5.90719 

3.37718 

1.85719 

7.58719 

8.8071 

2.63jri8 

1.20719 

4.99719 

3.38718 

1.55719 

6.41719 

8.9071 

2.64^18 

9.98718 

4.22719 

3.38718 

1.28719 

5.42719 

9.0071 

2.60^18 

8.25718 

3.56719 

3.34718 

1.06719 

4.57719 

9.1071 

2.S2IT18 

6.80718 

3.00719 

3.23718 

8.74718 

3.86719 

9.2071 

2.42iT18 

5.60718 

2.53719 

3.11718 

7.19718 

3.25719 

9.3071 

2.31^18 

4.61718 

2.13719 

2.97lri8 

5.92718 

2.74719 

9.4071 

2.18/T18 

3.78718 

1.79719 

2.80718 

4.86718 

2.30719 

9.5071 

2.02/?18 

3.11718 

1.51719 

. 2.60718 

3.99718 

1.94719 

9.6071 

).86/ri8 

2.55718 

1.27719 

2.39718 

3.27718 

1.63719 

9.7071 

1.70J718 

2.09718 

1.07719 

2.19718 

2.69718 

1.37719 

9.8071 

1.56F18 

1.72718 

9.01718 

2.00718 

2.22718 

1.16719 

9.9071 

1.42F1'! 

1.42718 

7.62718 

1.T2718 

1.83718 

9.79718 

1.0072 

1.20K18 

1.18718 

6.45718 

1.65718 

1.51718 

8.29718 

1.0172 

1.17^18 

9.78717 

5.48718 

1.50718 

1.26718 

7.04718 

1.0272 

1.06^18 1 

8.14717 

4.66718 

1.36718 

1.05718 

5.99718 

1.0372 

9.63K17 

6.80717 

3.98718 

1.24718 

8.73717 

5.11718 

1.0472 

8.76^17 

5.70717 

3.41718 

1.1371R 

7.32717 

4.38718 

1.0572 

7.99^17 

4.80717 

2.93718 

1.03718 

6.16717 

3.76718 

1.0672 

7*30f?17 

4.05717 

2.52718 

9.37717 

5.21717 

3.24718 

1.0772 

6,68^17 

3.44717 

2.18718 

8.58717 

4.41717 

2.80718 

1.0872 

6,13fl7 

2.93717 

1.90718 

7.C7717 

3.76717 

2.43718 

1.0972 

5.63K17 

2.50717 

1.65718 

7.23717 

3.21717 

2.12718 

1.1072 

5tlOFl7 

2.15717 

1.45718 

6.66717 

2.76717 

1.86718 
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Table A-4^ 


Alt 

— 

X = 60“ 



X=75» 







6.00^1 

4.22K18 

2.23921 

8.31921 

; . 69918 

4.25ni 

l*59ff22 

e.ioFi 

4-3in8 

1.94921 

7.24K21 

7.85P18 

3.71F21 

1.38ff22 

6.20ffl 

4.38918 

1.69921 

6 . 3092 1 

8.00918 

3.23F21 

1 . 20 n 2 

6.30ffl 

4.4S918 

1.47K1 

5.48921 

8.13918 

2.80921 

1.05K2 

6.40^1 

4.51^18 

1.27K1 

4.75E21 

8.25918 

2.43921 

9.07921 

e.soffi 

4.S7918 

1*10K1 

4.10921 

8.37918 

2.10921 

7.84921 

e.soffi 

4.82918 

9.47920 

3.54921 

8.47Jri8 

1.81K1 

6.77JP21 

6.70V1 

4.66ffl8 

8.1SF20 

3.05m 

8.57918 

1.56m 

5.83921 

6.80V1 

4.71918 

7.00920 

2.62F21 

8.66918 

i.34ni 

5-01921 

6.90V1 

4.7S918 

,S. 99920 

2-24921 

8.74918 

1.15921 

4.29ni 

7.00F1 

h.lBElB 

S. 11920 

1.91921 

8.82918 

9.78920 

3.66m 

7.10ifl 

4.82918 

4.35920 

1.63921 

8.89P18 

8.32920 

3.12921 

7.20FI 

4.8S918 

3.70920 

1.39ni 

8.96918 

7.osno 

2 . 66 m 

7.30fft 

4.389x8 

3.14920 

1.18921 

9.02918 

6.01920 

2.26^21 

7.40ffl 

4.90918 

2.66F20 

1.00921 

9.09918 

5.10920 

1.92m 

7. son 

4.93918 

2.'^5JT20 

8. 50920 

' 9.14918 

4.3in0 

1.63921 

7.60F1 

4.9S918 

1.90920 

7.2ono 

9.20918 

3.65K0 

1.38921 

7 •70^1 

4.97918 

1.61920 

6.10920 

9*25F1B 

3.08920 

1.17921 

7.B0F1 

4.39918 

i.seno 

5.16920 

9.30918 

2 . 6 ono 

9.8sno 

7 ,90^1 

S. 01918 

1.1S920 

4.37920 

9.35918 

2.19920 

8.37^20 

s.oon 

S. 03918 

9.66919 

3.70920 

9.40918 

1.85920 

7.08f20 

8. ion 

S.0S918 

8.13919 

3.13n0 

9.44918 

i.seno 

6 « 00^20 

8 . 20 n 

S. 08918 

6.84919 

2.65920 

9.49918 

i.3ino 

5.08^20 

8. son 

S. 08918 

5.76919 

2.25920 

9.53918 

l.^OWO 

4.30920 

e.4on 

S. 10918 

4.83919 

1.90920 

9,57F18 

9.26919 

3.64^20 

8. son 

S.llKlS 

4.0S919 

1.61920 

9.61918 

7.77919 

3.08f20 

8. eon 

S. 12918 

3.40919 

1.37no 

9.65918 

6.50919 

2.62n0 

8.70n 

S. 14918 

2.84919 

i.ieno 

9.69918 

5.43919 

2.77E20 

8. son 

S.1S918 

2.36919 

9.81919 

9.73918 

4.53919 

1.88120 

8. son j 

S. 18918 

1.96919 

8.29919 

9.76918 

3.76919 

1.58^20 

s.oon I 

S. 09918 

1.62919 

6.99919 

9.65918 

3.11919 

1.34^20 

9. ion 

4.93918 

1.34919 

5.90919 

9.36918 

2.56919 

l*iW20 

9.20n 

4.7S918 

1.10919 

4.97919 

9.01918 

2.11919 

8.SU19 

9. son 

4.S3918 

9.05918 

4.19919 

8.60918 

1.73919 

a.OlflO 

9.H0n 

4.27918 

7.44918 

3.52919 

8.11918 

1.42919 

e.74fl9 

9. son 

3.96918 

6.10918 

2.96919 

7.53518 

1.17919 

5.67F18 

9. eon 

3.64918 

5.01918 

2.49919 

6.92918 

9.58918 

4.7S»19 

9.70n 

3.34918 

4.11918 

2.10919 

6.34918 

7.8'»918 

4.01fl9 

9.80n 

3.0S918 

3.39918 

1.77919 

5.80918 

6.40918 

3.38fl9 

9.90n 

2.78918 

2.80918 

1.50919 

5.28918 

5.35918 

2.86fl9 

i.oon 

2.S2918 

2.31918 

1.27919 

4.79918 

4.43918 

2.42118 

i.oin 

2.29918 

1.92918 

1.08919 

4.34918 

3.67918 

2.06719 

i. 02 n 

2.08918 

1.60918 

9.16918 

3.94918 

3.06918 

1.75719 

1.03S2 

1.89918 

1.33918 

7.81918 ! 

3.58918 

2.55918 

1.49719 

i.04n 

1.72918 

1.12918 

6.69918 ! 

3.26918 

2.14918 

1.28719 

i.osn 

1.S7918 

9.42917 

5,75918 

2.97918 

1.80918 

1.10719 

i.oen 

1.43918 

7 96917 

4.95918 

2.71918 

1.52918 

9.45718 

1.07B7 

1.31918 

P.75917 

4.28918 

2.48918 

1.29918 

6.17718 

I.osn 

1.20918 

5.74917 

3.72918 

2.27918 

1.10918 

7.08718 

I.osn 

1.10918 

4.91917 

3.2A9te 

2.09918 

9.37917 

6.17718 

i.ion 

1.02918 

4.22917 j 

2.8S918 

1.92918 

8.05917 

5.40718 
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APPENDIX A-5 

Q-Function CompamenAs and r.m«s Error 


Calculated values of the total Q-fuiK:tion together with the important com- 
ponents are listed in tables A-5 for several values of solar zenith an(do> The 
final column labeled AQ/Q is our estimate of the relative error on Q(z, x ). To 
obtain this estimate we consider each compcment q. to have an independent rela- 
tive error given by 



/3FjV /3N-V 

-j— ) +[ !) +( li) +(1+(T..S.) 

vj \Ni/ Ui/ 



To evaluate this expression we assume somewhat arbitrarily the following 
uncertainties, 

flux = 25 

cross section ~ — ^ = .10 

* J * J y 

atmospheric density = .20 

nitric oxide 3NO/NO =1.0 

For the following components of Q we have used the uncertainties 


q(1025-991A) 

3q/q = 1.0 

q(977A) 

Sq/q = 1.0 

q(o*) 

9q/q = 2.0 

q(background) 

3q/q= 1.0 


The relative error on Q is then obtained from 



X 



2 
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TABLE A-5.1 


ra c. 

c > 


<i: < 

fc: 


5 ^ 
C 

C( r< fV. «« «« a4 

COOOOOOOOOOOOOOOOOOOCcjOOOO**r'^»-»*^*<**<.CNfJC»r»nr;«^r'»^-r: t- (*- nr; r> r> 

1 t 1 t 1 II OOOQCCOO'^OOOOni-r. rr'OOOC.CCOOOOOOOCOCCCCOCC OO O O 
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APPENDIX A-6 

Equatorial Reference Electron Density Profile 


In principal, the electron density distribution of the lower ionosphere from 
60 km to 105 km can be accurately deduced from ground-based radio wave ab- 
sorption and virtual height measurements, provided that the measurements are 
performed over a sufficiently wide range of freq'xencies, perhaps from about 
0.25 MHz to 3.5 MHz. In practice, it has provec to be difficult to conduct- experi- 
ments at frequencies less than 1.5 MHz partly because these frequencies fall 
within the broadcast band and also because the noon absorption is hig^ at these 
frequencies, so that the reflected wave cannot be detected by conventional methods. 
The frequencies generally employed for D- and E-region absorption studies 
have therefore been limited to the range 1.6-3.5 MHz. 

Electron density profiles have been derived from multi -frequency data 
acquired within this limiter’ frequency range (Beynon and Rangaswamy, 1969); 
but there seems to be some uncertainty in the precise form of these profiles at 
altitudes below about 80 km. This is because the total absorption and its varia- 
tion with frequency in the range 1.6-3. 5 MHz is rather insensitive to the exact 
shape of the profile at these low altitudes. Here the absorption is entirely non- 
deviative and is a relatively amall fraction of the total absorption. 

Ground-based methods which yield more accurate profiles for altitudes below 
80 km are the partial-reflecti-m technique, wave-interaction experiments and 
propagation measurements at low- and very-low frequencies. But none of 
these methods has hitherto been exploited at equatorial latitudes. The few 
profiles tnat are available for these altitudes at the equator are tlmse obtained 
from rocket-borne measurements at Thuirba, India (Kane 1969, 1972, Somayajulu 
et al. 1971, Aikin et al. 1972). 

In the present work we construct a reference electron density profile for a 
solar zenith angle X = 10**, for the conditions of an undisturbed sun (1-8A flux 
less than 0.5 x io~^ erg cm~^ sec'^ ) at solar cycle maximiim, by combining the 
results of a noontime rocket measurement at Thumbs, India with the results of 
multifrequency absorption and virtual heigdit’ measurements made at Colombo, 
Ceylon. For the lowest segment of the profile from 60 km to 82 km we adopt the 
electron density distribution derived by Kane (1969) from rocket-bome absorption 
measurements obtained when the solar zenith ai^le was 14**. This is the smallest 
zenith angle for which rocket measurements are available. The eleotrm density 
values computed by Kane have been corrected so as to be consistent with the 
current best estimates of collision frequency. 
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The ui^rmost portion of the profile depicts the E-li^er maximum. The 
corresponding electron densit 3 ' Is calculated from the critical frequency foE of 
the E-layer. Empirical formulas which give foE as a function of solar activity 
and zenith angle have been deduced from several decades of ionosonde data. 
The formulas used in the present work are listed below: 

(a) Minms (1964) 

For noon foE, 


where 


(foE)^ sec Xo 1.55 $ - 1.0 x lO'^^* + 14 


$ is the 10.7 cm flux in 10'^* Wm"* 

(b) Lyon (1965) 

Annual mean foE at Ibadan (0T‘22'N 03^ 58'E) around noon is given by 


foE= 3.47(1 + 0.0014R,) 


where R , is the sunspot number. Since the noon foE varies with season in pro- 
portion to cos^'^^ X this formula reduces to 

foE = 3.52 cos® ” x(l + 0.0014 R,) 


(c) Eyfrig and Rawer (1969) 


foE = 


3.44{1 - 0.1(l-co 



($ - 69) 



Cosi/3 y 


where 6 is the geographic latitude and is the 10.7 cm solar flux. 

For the solar conditions appropriate to our reference profile, that is 
v= lO**, R^ = 85 and $ = 135, the mean value of foE given by these formulas 
is 3.89 MHz corresponding to an elootron density of 1.9 x 10’ onf’ . 


The altitude of the E-layer peak Is calculated from the expression (Piggott 
and Thraiie 1966) 


h|||E = 105 + 7.5 ln(sec X) km. 
For y = 10“, 1\,E = 105.1 ^ 105 km. 
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Noontime rocket-bnrne probe measurements conducted at equatorial latitudes 
(Aikii'. and Blumle 1968, Satya Prakash et al. 1971) ha' ' indicated that the E-region 
has a broad maximum of Ionization which extends several kilometers above the 
peak altitude deduced from lonosonde records. For the reference profile we 
have therefore adopted a constant electron density of 1.9 x lo^ cm~^ for the 
altitude range 105-110 km. 

To determine the electron dejslty distribution in the intervening region from 
82 km to 105 km we first postula .e an arbitrary distribution and then adjust it by 
trial and error until the calculated values of absorption and virtual height of 
vertically incident radiowaves of different frequencies are in goo(' preement 
with the ground-basM measurements made at Colombo. The experimental data 
used for this purpose are those obtained with radio waves reflected from ffie 
E-layer, on the five frequencies 1.S3, 2.0, 2.2, 2.6 and 3.2 MHz. The data on 
1.33 MHz had been obtained by means of a special technique which enabled the 
detection of weak echoes (Gnanaiingam 1954). 

It is well known that the irirtual hci0it of reflection is particularly seimitive 
to the election density gradient near tlw reflection level. So, our initial trials 
aie directed towards finding a suitable profile for the bottomside of the E-layer. 

As a first step it is convenient to assume an electron density distribution that 
could be represented iiy one of the several formulas which have been postulated 
for the profiles of ionospheric layers (Budden 1961). The virtual helght onald 
then be readily computed from well-established analytical expressions, llie 
profile selected by us as a first trial is the 'sech ' profile given by the 
expression: - 


/h - h_ 

= 1.9 X 10* sech^ ( ■ 


where the symbols have their usurl meaning. When the earth's magnetic field 
and electron collisions may be neglected the virtual helf^t resulting from this 
profile is given by (Budden 1961, p. 157) 



h'(f) = H log 



This relation is also valid for an ordinary wave jted vertically near the 
magnetic equator. For H = 8 km, the virtual helgw. ^ computed from this 
expression are listed together with the measured valued; 
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Frequency (MHz); 

1.33 

2.0 

2.2 

2.6 

3.2 

True height (km): 

91.1 

94.7 

95.6 

97.3 

99.8 

Computed h' (km); 

96.8 

100.8 

101.9 

104.1 

107.8 

Measuv^ h' (km); 

96.5 

100.4 

101.3 

102.4 

104.6 


The discrepancies between the computed and me$isur~d virtual tieights, 
particularly at the two higher frequencies, indicate that the electron densities 
have to be increased slightly at altitudes between 95 km and 105 km. The profile 
is modified accordingly. It is also found that the lower end of the 'a^ch^' profile 
has to be rounded off so that it merges smoothly with the Kane profile at 82 km. 

We now have a composite profile from tO km to 110 km. The absorption «:nd 
virtu.-\ height resulting from this profile !<re next calculated by means of the 
y moralized magnetionic theory using the lay theory sq^nroximation. In calculating 
these quantities, the electron density and collision frequency values specified 
at 1 km intervals are logarithmically interpolated bv tlie computer for altitude 
laminations of 100 meters. Since the complex refractive index chuiges vari 
rapidly Immediately below the T*aflect; ' leve' , the last 1 km of the ray trajectory 
is subdivided into smaller laminations of 10 meter thickness. Tl» reflectica 
level h^ is defined for the ordinary mode as the altitude at which the plasma 
frequency equals the frequency of the exploring radiowave. 

To the absorption calculated by ray theory we add a phase integral correc- 
tion which is computed from the approximate analytic expression derived by 
Thorps (1971). Near the magnetic equator, this correction is relatively small. 
Thorpe has checked h - results against the phase integral correction computed 
rigorously, and has found that his accuracy is better than 5-10% provided that 
the angle between the wave normal and file goomagnetic field is greater fiisn 10° 
which, of course, is the case near the magnetic dip equator. At our ex;,loring 
frequencies, this combination of a ray theory caiailation plus a j^Aase integral 
co»*r.' ion gives a good approximation to a full wave solution. We estimate that 
the overall accuracy of the absorption calculated by this method is better than 
xl.O db. 


To calculate the virtual height we note that, for an ordinary wave propagat rg 
vertically near the magnetic dip equator, the phase refractive index n and the 
group refractive index n' ate related, to a good approximation, by the expression 
nn' - i, provided that the effect of colllrions may be neglected. Hence, from ray 
theory, the virtual height can be computed, with sufficient accuracy for our pur- 
pose, from the expression^ 


h' 


Jo 


dz 
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The advantage of using this expression is that n can be calculated more readily 
than n'. To the value of virtual hei^t calculated by ray theory we must add a 
I^ase integral correction so as to give a result that is closely equivalent to a 
full-wave solution. Instead however, for ease of computation, we adopt an alter- 
native procedure (Titheridge, 1967) in wuidi n(z) in the above int^ral is eval- 
uated with the collisi<» frequency set equal to aero. The accuracy of the virtual 
height resulting from this procedure is estimated to be better than ^0.5 km. 

The values <tf idteorption and virtual computed for the composite pro- 
file, for the five frequencies in the range 1.33-3.2 ^z, are now compared with 
the measurements. The profile in the altitude range 82-105 km is tiien sqitfdbly 
adjusted in file U^t oi any discrepancies that are noted and the above procedure 
is repeated. Th<s process of iteratiaii is continued until there is good agree- 
ment between the calculated and measured quantities. 

It is worth mentioning fiiat the values of idbsorpticm and virtual hei^ 
employed in the construction the reference profile are the mean values 
the noon measurements obtaiimd oa 72 days when the sun was umfisturbed, that 
is l-sA solar X-ray flux less than 5 x lO^ ergs cm~^ sec~*, in the equinoctial 
monfiis of 1968, 1969 and 1970 around solar cycle maximum. On 1.33 MHz, how- 
ever, measurements were only available for 16 days. 

The limits of uncertainty given in Table m for the meapured Misorpti<m and 
virtual height need some conunmat. The figures for absorptiqi)., Qiat is ±6 ^ for 
1.33 MHz and ±2 db for 2.0, 2.2, 2.6 and 3Jt MHz, are estimatra of file uncer- 
tainties in the calibrati<m cmistants, and hence in the absolute magnitudes the 
mean absorptimi. These estimates are based mi the number of ind^iendent 
measurements that were obtained for the calibratimi constant on each frequent, 
and (m the uncertainty (about ±1 db) in the ground reflectimi coefiicieat. 

The values (rf mean absorption are also subject to random errors arising 
partly from experimental errors and partly from the iidiereifi variMdUty of the 
ionosidiere. The precision of the individual absorption measurements is esti- 
mated to be better than ±2 db, and the standard deviation of the measuremmits 
in the 72-day sample is found to be idxxit 2.5 db. Hence, the standa.'^l error of 
file mean absorptimi is approximately {(2^ + 2.5^)/71}^'^^ = 0.38 db. This error 
is relatively small and is not therefore included in file lable. 

In the case of ihe virtual heigdit measurements, the calibratimi accuracy of 
the measurin'' system is better than ±0.1 km. Errors due to ixilse rise-time 
effects are es, nated to be less than 1 km. The principal source of uncertainty 
in the measurements is the distortion of the echo pulse caused by the 'roughness* 
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of the reflecting layer. This is a prominent feature of E-region echoes received 
near the magnetic dip equator and is a consequence of the irregularities associ- 
ated with the electrojet. The leading edge of the echo pulse is ill-defined, so that 
the virtual height cannot be determined with great accuracy. The reading accu- 
racy of each observation is no be^r than ±2 km. The iq^rent variation in the 
virtual heights measured on the undisturbed days is probably due to this cause: 
and tile standard deviation of these measurements may be regarded as indicative 
of the degree of rou^iness near the reflection level. 

In Table m, the uncertainties in tiie measured virtual height are represented 
by their standard deviations. From a statistical point of view the 95% confidence 
limits of the mean virtual heights are about one fourtii the standard deviations. 
However, considering all possible sources of error, the actual limits of uncer- 
tainty are prdiiably about :^1.5 km. 
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APPENDIX A-7 
Calculated (z, y) Profiles 

From a knoadedge of bow the Q-function changes with solar zenith angle 
(Appendix A-5) together with the noontime initial condition on the electron density 
profile (i.e. (ref). Appendix A-6) we can calculate the decay of the electi'on 
density profile as the solar zenith anj^e proceeds toward sunset. This we do 
by considering the following photochemical equati<m for the rate of chai^ of 
positive ion density. 


dN^dt = Q - - a.N^N. 


( 1 ) 


adiere , a. are the coefficients of electron-ion and ion-ion recombination 
respectively. Since 

N* = \ + N. = (1 + N, (2) 

dN^/dt (1 + \) dN/dt + N^d\/dt W 

Defining s (1 + \) (a^+ Ka.)» we then have 

(1 + K) dN/dt + N,d\/dt = 0 - 

If we now assume dhat \ ^ N_/N^ is time (zenith angle) independent, we 
have dX/dt = 0 and equation (4) becomes 


(1 +\)dN/dt=Q-a„N2 


( 5 ) 


^cre , is now time independent. This quantity can be evaluated by n<^ng 
that at noontime (i.e. X = 10**) dN /dt = 0, and therefore 


= Q(10«)/N^(ref) 
For the special case \ < < 1, equation (5) becomes 


( 6 ) 


dN/dt = Q-a„N’ 


( 7 ) 
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n^ch can be solyed numerically. When this is done we find that the dl^ /dt term 
is completely negligible for v < 75*. The decay of (z) with solar zenith angle 
can then be expressed as 


Ne(2- v) = N,(ref) 


/ 


Q(*» y) 
Q(*. 10*) 


( 8 ) 


The condition 1 breaks down at some altitude z 80 km. However 
since we have already deterir rted fimt the dN /dt term is ae^igihle for 1, 
it immediately follows from equation (5) fiiat w^/dt is even smaller when 
\ > 1. Thus within our assumption dX/dt = 0, equation (8) remains a good 
approximaticni at the lower altitu^. 

Profiles of N,(z,y) calculated from equation are shown in Fig. A-7.1,2,3. 
The upper and lower limits on N,(z,y ) shown in these figures result from the 
r.m.s. error A Q/Q given in Appendix A-5. Included in each figure is the <z) 
profile required above 95 km to yield agieement with the ground-based measure- 
ments of radiowave absoiption and virtual hei^. It is seen that the required 
profiles are well outside the error limits on our calculated profiles. 
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Figure A-7. 1. At x ^ required electron density profile compared with error limits on our calculated profile. 










REQUIRED PROFILE 
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At X » 75^ r«quir«d el#ction dMsily profiU compor#d with •rror limits on our colculoted proti 







APPENDIX A-8 
Equatorial N, (z,y) Model 

Tabulated values of our model for the afternoon decay of the electron density 
profile (shown in Fig. 8, Section 10) are given here in Table A-8. It should be 
noted that this equatorial model is for equinoctial months uxxler quiet sun con- 
ditions at solar cycle maximum. It should also be noted that at toe lower-most 
altitudes our model underestimates the decay of the electron density profile since 
we have assumed that the efiective recombination coefficient does ncrt increase 
with solar zenith angle (i.e.d\/dt = 0. 


PRECEDING PAG£ BLANK NOT FILMED 
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TABLE A-8 


Alt 

(km) 

Ne(10‘‘) 

(cm"®) 

N.(40‘*) 

(cm"®) 

N.(60') 

(cm-®) 

N.(75*) 

(cm"®) 

6.10F1 

7.50E1 

4.53E1 

4.26E1 

4.26E1 

6.20F1 

8.50E1 

4.00E1 

3.08F1 

3.08F1 

6.30F1 

9. 50F1 

4.12E1 

2 .24F1 

2.24fl 

6.40J71 

1.10E2 

4.98EI 

1.92E1 

1.82fl 

6. SOffl 

1.30E2 

6.41F1 

2.01F1 

1.70i?l 

6.60F1 

1. 55E2 

8.37E1 

2.54F1 

1.79F1 

6.70F1 

1.90E2 

1.11E2 

3.63E1 

2.07F\ 

6.80F1 

2.30E2 

1 . 4 5E2 

5,32F1 

2.50F1 

6.90E1 

2.90E2 

1.96F2 

8.09F1 

3.26E1 

7-OOFl 

3.50E2 

2.50F2 

1.16E2 

4.18E1 

7.l0ffl 

4.20E2 

3.15F2 

1.63F2 

5.53E1 

7.2GF1 

5.00E2 

3.92F2 

2 .23F2 

7.58E1 

7.30F1 

5.70E2 

4.63F2 

2.85F2 

1.03E2 

7.40ffl 

6.30E2 

5.27E2 

3,49F2 

1.38E2 

7.50B’l 

6,80E2 

5.84E2 

4.10F2 

1. 80E2 

7.60F1 

7.40E2 

6.49E2 

4. 80E2 

2.33E2 

7.70F1 

7.80E2 

6.97F2 

5.38F2 

2.88E2 

7.80E1 

8.40E2 

7.61E2 

6.10E2 

3.S5E2 

7.90F1 

9.80E2 

8. 07F2 

6.66F2 

4.19E2 

B.OOffl 

9.10E2 

8.42E2 

7.13E2 

4.78E2 

8.10E1 

9. 30E2 

8.67F2 

7.49E2 

5.31E2 

8.20E1 

9.50E2 

8. 89E2 

7.81F2 

5.80E2 

6. 30F1 

9.85E2 

9.23E2 

8.22F2 

6.36E2 

8.40E1 

1.07E3 

1.00F3 

9.00F2 

7.21E2 

8.50F1 

1.30E3 

1.21E3 

1.10F3 

9.06E2 

8.60E1 

1.85E3 

1.70E3 

1.55E3 

1.31E3 

8.70E1 

2.75E3 

2.48E3 

2.25E3 

1.95E3 

8.80F1 

4. 0OE3 

3.48F3 

3.14F3 

2 .27E3 

8.9 0/?l 

5. 80E3 

4.86E3 

4.32F3 

3.88E3 

9.00ffl 

8.20E3 

6.61E3 I 

5.71E3 

5.18E3 

9.10J71 

1.12E4 

8.62E3 

7.11F3 

6.47E3 

9.20E1 

i.55E4 

1.1 3F4 

8.78F3 

7.93E3 

9.30ffl 

2.10E4 

1.44E4 

1 . 04F4 

9.17E3 

9.40ffl 

2.70E4 

1.76E4 

1.1 5E4 

9.78E3 

9.50/?l 

3.55F4 

2.25E4 

1.30F4 ! 

1.06E4 

9.60ffl 

4.75E4 

2.60E4 

1.42F4 

1.08E4 

9.70E1 ! 

6 . 2 0E4 

3.10E4 

1.55F4 1 

1.10E4 

9.80P1 

8.00E4 

3.90E4 

1.80F4 

1.12E4 

9.90E1 

1.00E5 

5.00E4 

2.15E4 1 

1.15E4 

1.00F2 

1.22E5 

6.60F4 

2.70F4 

1.21E4 

1.01F2 1 

1.44F5 

8.60E4 

3.80F4 i 

1.32E4 

1.02F2 

1.62E5 

1.07F5 

5.00F4 

1.50E4 

1.03F2 

1.76E5 

1.25F5 

G.20F4 ! 

1.80E4 

1.04F2 

1.85E5 

1.40E5 

7.20F4 

2.50E4 

1.05E2 

1.90E5 i 

1.51E5 

8.50F4 

3.90E4 

1.06E2 

1.90E5 

1.S0E5 

9.70F4 

5.60E4> 

1.07F2 

1.90E5 

1.65F5 

1.08F5 

6.50E4 

1.08F2 

1.90E5 

1.65F5 

1.15F5 

7 . . 3E4 

1.09ff2 

1.90E5 

1.65E5 

1.20E5 

7.30E4 
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